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Abstract 

NASA is planning missions to the vicinity of the Sun-Earth L 2 point, 
some involving a distributed system of telescope spacecraft , configured in 
a plane about a hub. Several sets of differential equations are written for 
the formation flight- of such telescopes relative to the hub, with varying 
levels of fidelity. Effects are cast as additive perturbations to the circular 
restricted three-body problem, expanded in terms of the system distances, 
to an accuracy of 10-20 m. These include Earth's orbital eccentricity, lunar 
motion, solar radiation pressure, and small thrusting forces. Simulations 
validating the expanded differential equations are presented. 

INTRODUCTION 

NASA's Goddard Space Flight Center is planning a series of missions in the vicinity of 
the Sun-Earth L2 libration point. Some of these projects will involve a distributed space 
system of telescope spacecraft- acting together as a single telescope for high-resolution. The 
individual telescopes will be configured in a plane, surrounding a hub. where the telescope 
plane can be aimed toward various astronomical targets of interest. Nominal missions locate 
the hub in a halo orbit, at distances of hundreds of thousands of kilometers from the L2 
point. The individual telescopes are to be placed at less than one kilometer from the hub. 

The present research continues the earlier work of Segerman and Zedd [l]. In that 
paper, the circular restricted three-body problem was used as a basis for the derivation of 
differential equations which describe the motion of a typical telescope spacecraft relative 
to a nearby hub. Following the approach of Richardson [2], a modified Lindstedt-Poincare 
method was used in order to develop an uncontrolled periodic solution for this relative 
motion, which is a requirement for maintaining the telescope formation. The solution 
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included linear effects of the hub motion about Lo: quadratic hub motion effects were used 
to develop relationships between the telescope frequency and amplitude. In the course of 
that analysis, a halo-type orbit of the telescope about the hub was chosen to provide periodic 
motion in the aperture plane. Comparison was made between the resulting relative motion 
solution and a baseline numerical solution of the full circular restricted equations of motion. 

This research furthers the earlier work by adding the increased fidelity associated with 
perturbations to the circular restricted three-body problem. Each of the perturbing effects 
is cast as additive perturbations to the circular restricted equations of motion. For the 
purpose of this analysis, these perturbations include the ellipticity of Earth's orbit about 
Sun, the orbiting Moon, and solar radiation pressure. 

In considering the effects of the elliptical restricted problem, it is noted that the various 
relative motion vectors within the three-body system are ordinarily functional!)' dependent 
upon Earth's orbital eccentricity. To avoid this dependence, the vectors are denned relative 
to reference positions of Sun. Earth, and Li along the line of syzvgy. Then, the relative 
equations of motion for an individual telescope are expanded in terms of the distance of the 
hub from Li as well as the distance of the telescope from the hub. These expansions are 
written so as to include all terms which contribute at the 10-20 meter level or above. Explicit 
linear contributions of Earth's orbital eccentricity are retained; higher order contributions 
are excluded, as being below the prescribed level of fidelity. 

For the inclusion of the lunar contribution, expansions of the vector forces are formed, 
here involving the small-scale distance of a telescope from the hub, as compared to the larger 
scale of distances of the lunar motion within the Sun-Earth system. Again, the expansions 
are truncated in a fashion consistent with the previously stated fidelity level. 

Solar radiation contributions are included as well. These terms include the effects of 
possible eclipses resulting from various geometric configurations of the system. This effect is 
included both for completeness and with the awareness that it may be possible to actively 
adjust the spacecraft attitude in order to use the solar radiation pressure as part of an 
orbital control scheme. 

Finally, the equations are written to include the presence of small body-fixed thrusting 
forces, which are assumed to be used as part of the orbit control. 

Simulations were conducted, comparing numerical integration. of the various resulting 
sets of differential equations with numerical integration of the full, unexpanded differential 
equations. These simulations have been used to verify the validity of the expanded equa- 
tions. Some of the results of these comparisons are presented. Studies of the solution’s 
sensitivity to errors in hub position knowledge are also discussed. 


EQUATIONS OF MOTION 

Circular Restricted Three-Body Problem 

In the preliminary phase of this research [1], the general second order differential equations 
of motion were constructed for an object near the Sun-Earth Li libration point, using the 
force model of the classical circular restricted three-body problem. In this model, Earth is 
treated as being in a circular orbit, about the sun, the spacecraft mass is considered to be 
negligible as compared to the two primaries, and only point-mass gravitational forces are 
considered. 

For this system, depicted in Figure 1, the differential equations of motion for an object 
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Figure 1 Coordinate Axis Definition 


(object i) near the Sun-Earth Lo are given by 


where 


r 7 — 
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x — rrr £ x« 


r ? = vector from Lo to object ?’ 
yi = solar Kepler ian constant 
fj -2 = terrestrial Keplerian constant (Earth — Moon) 
pn = distance from Sun to object i 
P2i = distance from Earth-Moon barycenter to object i 
x e = distance from system barycenter to Lo 
Di = distance from system barycenter to Sun 
Do — distance from system barycenter to Earth-Moon ban* center 
x = unit vector parallel to Sun-Earth line of svzygy, 
pointing in Sun-to-Earth direction 
n = terrestrial mean motion about Sim (assumed constant). 


The coordinate frame of Figure 1 is a rotating reference frame with origin O at the 
system barycenter. The maxis points along the Sun-Earth line of svzygy ana the z-axis 
is parallel to the system angular momentum: the y-axis completes the dextral coordinate 
system. 

Let r/j and r t denote, respectively, the vector from L 2 to the hub and to a telescope. 
Therefore, if r is the vector from the hub to the telescope, the differential equation of motion 
for the telescope relative to the hub is 


r = -u 1 
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where now. in general, the subscripts h and t refer to the hub and telescope. 

Elliptical Restricted Three-Body Problem 

The extension to Equation (1) is now developed for the case of the elliptical restricted three- 
body problem. First, it is necessary to locate the point which is analogous to the libration 
point L -2 in the circular restricted problem. As would be expected, such a point exists; as 
the Sun-Earth distance varies, the location of the point oscillates along the x-axis. 

Elliptical Problem Libration Point Analog. Say that there is an elliptical analog to 
the L 2 point and that its position relative to the (assumed inertial) system barycenter is 
given b} r 

R e = x e x 4- y e y - Z e Z. 

Letting / refer to the true anomaly of Earth’s orbit about the sun, the coordinate system 
has angular velocity 

Wc = /Z, 

which now is considered to be varying throughout the year. Differentiating R fi , 


’ Xe - he ’ 


’ Xe - he ~ 2fy e - j 2 X e ' 

Ve 4 fx c 

and R^ = 

Ve + he + 2/x e - / 2 y e 

L J 




where the column vector notation is used to indicate the xyz vector components. 

‘The Newtonian gravitational force per unit mass acting on an object at this point by 
Sun and Earth is given by 


where 


F e /m = 


ui(x e 4 D\) /j 2 {X' : -D 2 ) 

o o 
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rie = (x e -T £>i)x + y e y 4 z e z 


and r 2e = (x e - Dn)x - y e y 4 z e z. 


and D\ and Do now refer to the time-varying locations of Sun and Earth along the line of 
syzvgy. 

Clearly, as in the circular restricted case, the z equation is decoupled, and admits a 
solution z e = 0. If, as anticipated, the desired solution involves y e = 0, the x and y 
components of the force- acceleration equation become: 


x: x e - f~x e -f 7 


Ml 
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From the y- component equation. x e ~f is therefore constant. However. conservation of the 
Sun-Earth two-body angular momentum h per unit mass sores 

D~f = h, 


where D is the varying Sun-Earth distance D\ 4 £> 2 * Of course, unlike the case of the 
circular restricted problem, this distance varies throughout the year. Substitution for / in 
the y- component equation gives 



constant. 


Fui laiei ounvciiiciitc. Jelliie this; coiistaiu. In i-erius of the constant 7, such that 



Taking the positive root, this gives the definition of 7 as 

_ h Ze - £2 

^ ~ D * 

In the case of the Sun-Earth system, 7 is approximately 0.01007824; x t varies between 
1.486 x 10 s and 1.536 x 10 s km throughout, the course of the year. 


Relative Motion Equation Derivation. Let R-; denote the position vector from the 
system b arv center (point O) to spacecraft i. Referring again to Figure 1 . 

R t = a? e x 4* r*. 

Differentiating, 

FL = x t ± 4 fx e y 4 r ; 

R; = (x e - fx e )x - (2 fx e 4 fx e ) y - r x . 

The distance x e and its time derivatives may be written in terms of the varying Sun- 
Earth distance D and its derivatives, using the following definitions of the constants 7 and 
p, and the associated relationships: 



A X e - D 2 

1 ~ D 

- A M 2 
V 

= pD 


, 1 Xe + A 

' ’ D 

1-/9= — 

M 

D'i = (1 - p)ZX 

where 


M = Ml 4 M2- 


Then. 

x e = 7 D 4 - D 2 
= (7 4 1 - p)D, 

x e = (7 4 1 — p)£?. 

Xe = (7 -4 1- 
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This gives the acceleration vector R,: as 

R, = ( 7 + i - p)[D - pD)x - (7 - 1 - p)(2fD - fD.ry -f r>. (2, 

Similarly, two-body relationships may be used to write D and D in terms of / and /. 
Using the two-bodv equations of motion for the sun (could use Earth): 

Ri = -D\k = -pDx 
Rj = -pDx - pfDy 
Ri = p(f 2 D - D)x - p(fD 4- 2fD)y 

n 2 

= —x (two-body force) 


This gives the component equations 

x:p(f 2 D-D) = -gr 

y '• -p(fD + 2/0) = = o, 


which yield 

D - j 2 D — -db and -/ D - fD = 0. 

pD~ 

Substituting into Equation (2) : the acceleration becomes 
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Introducing the two-body forces, 
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Note that the explicit appearance of the time derivatives of / has now been removed. 
Next, the vectors p Vl and p 2r may be written as 


Pu = (z e 4- D^x-rq 
— (7 -f l)Dx 4 r, 


and 


Substituting into Equation (3). 
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In the rotating frame, rq and its time derivatives are given by 
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Vote that vectors p—, pv : - ?; nd iy ore inactions of tiie eccentricity of Earth's orbit, as is 
D. To avoid this dependence, redefine these vectors relative to reference positions of Son. 
Earth, and Li alone' the line of syzy£.v. Denoting' the reference distances by an over bar. 

R ? = x t x — r ? 

Ri = fx e y “ r* 

R< = -/ 2 f e x -f / x e v — ib. 

Again using the definitions of 7 and p. x e = (7 -f 1 — p)Zh 
Using this form of the acceleration. 


R, = (7 T 1 - p){— f 2 x -r fy)D — iy 
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where 

Piz — (7 T l)Px 4- r, and p 2l = 7 Dx ~ r*. 

Vectors r- t and iy take the same forms as before, now relative to the reference Lo location. 
For the relative motion, let 


r = r t - r h = R* - R,v 


Then, 


r = -pi 


rt 
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Note the similar form to Equation (1). 
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Lunar Gravitational Effects 

This section discusses the lunar contribution to the telescope motion relative to the hub. 
Terms corresponding to the gravitational force of the moon upon the hub and telescope 
spacecraft are included in the relative equations of motion (telescope relative to hub). The 
resulting contribution is then cast as an additive perturbation to the elliptical restricted 
three-body problem equations, such that when the lunar motion is ignored, the contribution 
reverts to the lunar mass placed at the Earth position. The moon is treated as a point mass. 

As shown m Figure 2, let denote the vector trom tire moon to spacecraft i, either 
the hub or a telescope, and let 73 denote the lunar Keplerian constant. The lunar force per 
unit mass on spacecraft i is then given by 


Pa i 

-M3 3- 

P 3 i 6 


( 6 ) 


Therefore, the contribution to the equations of motion of the telescope relative to the hub 


is 


_ ... f _ EMl ' 
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■where the subscripts h and i refer respectively to the hub and telescope. 
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Figure 2 Coordinate Axis Definition (including Moon) 


Solar Radiation Pressure Effects 


The pressure from solar radiation imparts a relatively tiny force on a telescope spacecraft. 
Depending upon spacecraft design and distance from the sun. the force can perturb both 
the spacecraft's attitude and orbit. The force can also be harnessed to beneficially propel 
the spacecraft. The model presented here provides the force on a spacecraft, accounting for 
the force reduction when the spacecraft orbits through the terrestrial shadow*. 

At a distance p\ t from the sun, the solar flux / (the irradiance) acting on the spacecraft 
is given by 


/ = 


L 

^ Pit 


where 

L - 3.842 x 10 2G watts 


is the solar luminosity (total emitted radiation). 

If A is the cross-sectional area of a spacecraft of mass m projected normal to the 
spacecraft- Sun line, then the solar radiation force F s per unit mass acting on the spacecraft 
is 


F, 


C r LA 1.0198 X 10 li C R A 


m p'i t 


N/unit, mass, 


A7imcpl t 

where c is the speed of light and 0 < Cr < 2 is the parameter characteristic of the reflectivity 
of the spacecraft surface facing the sun: 


Cr — 0 translucent, 

Cr — 1 perfectly absorbent, 
C R — 2 perfectly reflective. 


For a trajectory which passes through any portion of Earth’s shadow, the full disk of 
the sun will be partially obscured. In the vicinity of Ln. this will occur at distances normal 
to the line of syzygy of approximately 13,420 km or less. In such cases, the force expression 
above must be scaled by a “luminosity reduction factor" a which ranges from zero (total 
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eclipse; to unity "'full sunlight }. The appropriate expression for the 
then written: 

1 . 0 i 9 3 x 1 U * 1 CrA cr - . 
m \ I (7 -f 1 )DSc -t- r t j | 3 " 


F, = 


r- -i r, o 
U 1 3 O 


„ . r, : : .4 

a L pxr 

m plt 


■pi/ = 


force per unit mass is 
IjDX - Tt: . ( 3 ) 


Spacecraft Thruster Effects 

This section presents the incorporation of the effects of body-mounted thrusters in the 
equations of motion. 

First, the body-fixed acceleration components imparted by the thrusters are given as 


x b = Fx/'m, y b = Fy/m , z b = iT/m, 

where F y . F z are the components of the thrust F t > irust in the body- fixed frame: rn is the 
vehicle mass. 

Consider an arbitrary alignment of a body-fixed coordinate frame with respect to the 
rotating x. y, z frame. This body-fixed x> n yj , . z b coordinate frame has its origin at the 
spacecraft's mass center. 

The components of thrust expressed in the body frame must be transformed into the 
rotating reference frame in order to correctly incorporate these forces into the description 
of the motion. This transformation is expressed as 


X 


" X b “ 

v 

= T 

Vb 

Z b 


where T is the transformation matrix formed as a combination of Euler rotations. 

It is somewhat easier to visualize the individual rotations by considering the inverse 
rotation description, rotating instead from the rotating reference frame to that fixed in the 
spacecraft- By first assuming the coincident alignment of both the body frame and the 
rotating reference frame, the inverse transformation matrix T~ 1 is then a combination of 
Euler rotations through a set of Euler angles: 

1. first, rotate angle (ip -f /; about the z b axis, where / = f(t - to) 

2. second, rotate angle 9 about the new orientation of the y b axis 

3. third, rotate angle 6 about new orientation of the x b axis 


Combine the sequence of rotations as follows with a right-to-left ordering of the rotation 
matrices as 


'i 

1 

o 

o 


’c{9) 

0 

-swy 

0 

c(d») s{6) 


0 

1 

0 

_0 

-s(<p) c(6) 

i 

s{9) 

0 
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c{ib + f) s(ip + f) 0 
c{ij> + f) 0 
0 0 1 


where the functions c and s represent cosine and sine, respectively; the rotation matrix 
T-~ l refers to the required Euler rotation matrix about the f^-axis. Finally, the desired 
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transformation matrix T (from the body-fixed frame to rlie rotating frame) is expressed as 
the inverse of T” 1 : 

t = (T- l y l = r z (v - fjZjojTjo). 

These expressions premultiply the thrust force per unit mass acting upon the telescope 
spacecraft . 

Note that typical accelerations due to thrust and solar radiation pressure are several 
orders of magnitude smaller than those due to terrestrial and solar gravity. It is noteworthy 
that solar radiation pressure can be used for orbit control because the force from solar 
radiation pressure may be comparable to that from a thruster suite. 


Summary 

The force equations derived above are now combined. The resulting equation models the 
elliptical restricted three-body problem, incorporating lunar, solar radiation, and thrust 
perturbations. 

Combining the expressions of this section, given by Equations (4), (6). and (8), along 
with the thrust expressions above, the differential equations of motion of the telescope are 
given bv: 


Rt = (7 -t 1 - P)(-/ 2 * + fy)D - r t 
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where w f = v + / . The corresponding equations for the hub are 


R/i = (7 4- 1 - p)(-/ 2 x -r fy)D 4 r h 

f Mi 
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Combining the relative motion expressions of this section (Equations (5) arid (7)), along 
with the solar radiation pressure of Equation (8) and the earlier thrust expressions, the 
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In the MAXIM or similar missions, there may or may not be an actual hub spacecraft 
located at the aperture’s center. Regardless, it is necessary to treat the hub as a central 
reference point for locating the positions of the individual telescope spacecraft- However, 
unlike the gravitational forces present , the solar radiation pressure effect upon a spacecraft 
with actual mass and area is substantial as compared to that upon a "phantom" hub. 
Therefore, for purposes of simulating the solar radiation pressure, the hub “spacecraft” is 
treated as having the same physical characteristics as the telescope. In this manner, the 
hub is maintained as an adequate reference for the position of the telescope. 

An example is presented in order to permit comparison of the relative contribution of 
the terms to the telescope motion. The initial conditions are listed in Table 1; they were 
taken from the examples discussed in the earlier phase of the research. As mentioned in 
that paper, these values were selected so as to excite only the oscillatory linear modes. The 
relevant physical constants are presented in Table 4 below: the computed value of D may 
be found in Ta.ble 5. 

Figure 3 depicts the solution to numerical integration of four different force models 
selected from the summary equation above. The models were applied to both the hub and 
telescope, and the resulting state vectors were differenced in order to determine the position 
of the telescope relative to the hub. In each case, the same force model was applied to both 
the hub and telescope. Additionally, the same value of the reference distance D was used 
for all force models. 

The solutions are displat^ed relative to the circular restricted model solution, which is 
obtained by including the elliptical restricted model (ER3B) terms, but with Earth’s orbital 
eccentricity treated as being zero; this duplicates the model discussed in the earlier research. 
The other models depicted in Figure 3 represent the addition of ellipticity (ER3B), lunar, 
and solar radiation (SRP) effects over 20 days. It is noted that the elliptical contribution 
provides the dominant perturbation to the circular restricted solution. In this example, 
should the viewing of a scientific target be limited to no more than 10 days, the perturbation 
due to the elliptical contribution is less than 1 m. 

The effects of thrusters are not simulated here, due to the vast uncertainties of force 
(both magnitude and direction) and duration. 
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Table 1 

Example Initial Conditions 



hub 

(state rel. 
to L-o) 

telescope 
(state rel. 
to L 2 ) 

telescope 
(state rel. 
to hub) 

x(0) (km) 

-227,219.419 

-227,219.483 

-0.064780 

1 /( 0 ) (km) 

0.0 

0.0 

0.0 

z( 0) (km) 

-250,000.000 

-249,999.974 

0.026445 

i(0) (km/day) 

0.0 

0.0 

0.0 

y( 0) (km/day) 

25,625.039 

25,625.044 

0.004421 

i(0) (km/day) 

0.0 

0.0 

0.0 

mass m(kg) 

cn 

O 

O 

500 


^thrust vNj 

0 

0 


sun-facing area A (m 2 ) 

150 

150 




Figure 3 

Effects of Perturbations on Relative Distance — Full Equations 


EXPANDED EQUATIONS OF MOTION 


Circular Restricted Three-Body Problem 


In the previous research [1], Equation (1) was expanded through terms which are linear in 
the coordinates of r and no more than cubic in the coordinates of r/ } ,. 

The acceleration vector r may be written relative to a rotating coordinate system which 
rotates at the constant angular rate n about the z-axis normal to the ecliptic, and with the 
x direction as previously defined. This gives 

!~ x — 2m's — n 2 .r "1 

. 0 i 


r — 


4- 9. 


Znx — rr 


< 4 , 
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Elliptical Restricted Three-Body Problem 


O r.i' c; rjov fhp ^PctT'r-tPQ f lire^ - !' odv P C ■ iat ~4 '"wf^ 

The nglii tide ul Luis set of relative deceleration equal ions may 
research. Consider the effects of various contributions to the 
For ordering purposes, take 


o n as given i n E g u a r i o n , o : . 
be expanded as in the earlier 
magnitude ordering scheme. 


(r/,. = 600, 000 km) 

(r = 0.5 km). 

A rough estimate may be obtained of the contributions that- the various perturbing 
contributions make to the circular restricted problem solution. Say that a perturbing term 
may be treated as modifying the linear frequencies associated with the circular restricted 
problem, and consider the square of the perturbing frequency to be roughly the magnitude 
of the coefficient of r in the perturbing acceleration. (The linear periods in and out of the 
xy -plane are approximately 177.566 days and 184.002 days, respectively.) Then, after 90 
days, the effects of the terms containing various powers of r and r- n are given in Table 2, 
with effects included of roughly 20 in and larger. 


p = 3.04 x 10" 6 
7 = 1.01 x 10~ 2 
^ = 4.0 x 10~ 3 

V = 3.2 x]0" 9 
D 


Table 2 

Along-Ellipsoid Effect of Sun-Earth Perturbation on Solution (90 days) 


j perturbing 
1 term 

effect 

(m) 

j re 

! 16.7 

| rr h 

289.3 

| rr h 2 

118.5 

^ 3 
; TT h 

47.4 

rr h 4 

18.9 


Retaining terms which contribute to approximately 20 m, again through linear terms in 
the relative motion and cubic terms in the hub position, results in the expanded/truncated 
differential equations 


r = A 3 [— r 4 - 3xx] 

4- A 4 [3xr h 4- 3x h r 4- (Sr* • r - loxx/Jx] 

+ -4 3 [(3r A • r - 15xx h )r h -f \{r k 2 - ox h 2 )x - f{2x h r h • r - lxx h 2 -f xr h 2 )x] 
+ -46 [^rh{-xr h 2 - 2x /T ' r h - ’ixxk 2 ) 4- fr(7x ft . 3 - 3 x h r h 2 ) 

4- y(— r h 2 ri,. ■ r -r lxrf 2 Zh + 7x ;, 2 r ■ iy - 21xx/, 3 )x|. 


v;here the constants .4,, tabulated in Table 5. are given by 

4 . _ 44 _L_ tL 2 

(x e — Z>i) ! ' (x e — Do ) 1 


(9) 
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Next consider the left side of the differential equation: 


* - fy - 2 fy - f 2 x 

a - fx + 2jx - z 2 ?/ 


x — 2n c y - n c 2 x 
y - 2n c x - n c 2 y 


- fy ~ 2 (Z - n c )y - (J 2 - n ( r)x “ 
fx -r 2(/ - n c )x - (Z 2 - n c 2 )?/ 

0 


( 10 ) 


In this expanded form, the first vector term represents the acceleration which appears in 
the circular restricted problem (n c refers to the mean motion of the circular restricted 
Earth orbit). The second vector term gives the perturbation which is added by including 
the elliptic restricted effects. The perturbation term is to be expanded in terms of the 
eccentricity e of Earth's orbit (rs 0.017). In keeping with the earlier magnitude ordering, it 
is estimated that it is sufficient to retain only contributions which are at most linear in e, 
as er is approximately 8.5 m. 

Say that the circular restricted problem takes D as being the reference value D. Then, 
the corresponding mean motion is n c — ^ / \±f D 3 . Now. in the elliptic restricted problem, 
the mean motion is 


= la = .fjfjf 

V a 3 V D 3 a 3 



n c . 


If D is chosen to be the semi-major axis «, then n ~ n c . Alternatively, if D is taken as the 
mean Sun-Earth distance a(l -f e 2 /2), the linear eccentricity approximation still allows n to 
be approximated by n c . 

Using two-body relationships, and linearizing in e, the time-varying f is now approxi- 
mated by 

/ = ^2 ~ ~ n c(! + 2e COS /) 

or, in terms of the mean anomaly £ y 


f % n c ( 1 -r 2ecos£). 


Differentiating, 


f ^ -2 en c 2 sin A 


Using these results in the perturbing vector of Equation (10) gives 


' -JV - 2 (/ - n c )y - (/ 2 - n c 2 )x ' 


—2 eyn c 2 sin £ — 4 eyn c cos £ — 4 exn c 2 cos £ 

fx 4 - 2 (/ - n c )x - (/ 2 - n 2 )y 

= 

2 cxn c 2 sin £ ~ 4exn c cos £ - 4cyn c 2 cos £ 

o J 


0 


Lunar Gravitational Effects 


Consider the lunar force given above by Equation (7). It is desired that this contribution 
be expressed in terms of the position of the telescope relative to the hub. Let the vector r 


1:0 pUOJ I'lOii . 


Pzt ~~ Pzh “f r - 


( 11 ) 
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As in the earlier analysis of the effects of Sun and Earth, the co 
may be written as an expansion in terms of r and its comp one: 
binomial expansion development is followed. The square of the 


nrnbuT i u n oi E qu at ion ' / j 
its. Accordingly, a similar 
magnitude of p v is eiven 


bv 


P3t~ “ P3t * Pzt = PZk 




2 pzh ' r 
PZh~ 


Using this formulation to perform a binomial expansion of Equation (7) . using Equation ( 11 ) 
to substitute for p it) and expanding through linear terms in r gives 


F 


m 



-r 3/i 3 


Psh r ‘ P'Sh 
Poh 3 P3h 2 


(12) 


It is preferable to treat these terms as an additive perturbation to the equations of 
the elliptical restricted three-body problem. Assume that the baseline elliptical restricted 
problem contains an object with combined terrestrial and lunar mass, located at the mean 
Earth-Moon b ary center. Therefore, the expansion of Equation (7) should contain the lunar 
contribution to this combined mass along with terms representing the effect of the lunar 
motion about the bary center. Accordingly, Equation (12) may be rewritten as 




\P 2 h' J 


yPgA r ’ P2h \ 
P2h 3 Pll? ) 


- Pz* 




( PzhPzh 
\ PZk° 


P2hP2h \ 
P2h° ) ‘ 


(13) 


Because the vector p 2h used in the earlier analysis, it is convenient to write p 3h as 


Pzk — Pik T r em , 


where r em refers to the lunar position relative to the mean harycentcr. Using a Legendre 
polynomial expansion, the inverse of the magnitude of p 2h is given by 


1 


P3h 




P2h ' r c?T 

V P2h r em 


Expanding through terms linear in r cm , and substituting in Equation (13), 




P2 


, 3 


3p 2 /i r P2h 
P2h 5 


- 3p 3 \JP2h ' r em + t ■ p 2h r em A r • r emp2h 




Dr ‘ P2hP2h.P2h ' r em 


/P2fc 2 l /p2> 


’2 h 


Examination of Equation (14) indicates that the form of the first term in F m is identical 
to the corresponding terrestrial term derived in the earlier work. The only difference is that, 
here, the mass coefficient is ^.3 rather than p 2 . This term represents the effect of a lunar 
mass collocated with the terrestrial mass; the remainder of F m represents the effect of the 
lunar motion about Earth. The first term may be added to the earlier work simply by 
replacing p 2 with p 2 - p .3 in the relative equations of motion for the elliptical restricted 
three-bo dy problem. 
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For the second term of Equation (14). the Ea.rth-hub position vector p 2 h 13 expanded 
in a fashion similar to that of the earlier work. In the context ox the elliptical restricted 
problem using a reference location of Lo, 


p 2h = 7 Dx. -r v h 


where 7 and D are constants as previously defined. The square of the magnitude of p 2h is 
given by 


P2h 


fs2 


2 “ P2h ' P2h = 7-D 


rh 
7 D 


7 D 


where Xh is the 2 - component of r^. Once again, inverse powers of p 2 h are formed using 
binomial expansions. Expanding through linear terms in 77 , and substituting. F m becomes 


F 


m 


r _ 3p 2 fcf ' P2h 
P2h i P2h 5 

3/^3 |( 2 XXf> rn , T yyerri ~r ZZ eTn )yi -f- (jJX era ~r Xy ( . rn ) y 

+ [zTem + XZ em )z]/(yD) 4 
- / 2 3 [- lox em x h r + 3r/i • r em r - 15xr em 2: ft 

-i- 3r • r,,r em -f 3(r - v em ){-ox h x -f r„) + I0oxx cm x h x 
- lSarXcmr* - 15xr fc - r Bm x - 15x CTO r • r^xJ/OZ?) 0 . 



(15) 


Again, a rough estimate of the contributions that the various perturbing contributions 
make to the solution is presented. As was done for the inclusion of the ellipticity, the 
perturbing terms are treated as modifying the linear frequencies — that is, the coefficient 
of r in the perturbing acceleration. After 90 days, the effects of the dominant lunar terms, 
containing various powers of the relevant variables, are given in Table 3. 


Table 3 

Along-Ellipsoid Effect of Lunar Perturbation on Solution (90 days) 


perturbing 

effect 

term 

(m) 

rr em 

2.3 

9 

TV r JW „ 

0.6 

rri,r em J 

0.9 


Summary 

The force equations derived above are now combined, repeating the models for solar radia- 
tion pressure and thrust. The resulting equation models the expanded elliptical restricted 
three-body problem, incorporating lunar, solar radiation, and thrust perturbations. The 
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nous ox relative motion ror a relescoi 


spacecraft are given by: 


f = r — oxx - .4 ! 3xry- ( - 3v;,,r — f 3r/ t • r - 15 xxh.V 

- -4 5 ; Or;, • r - loxx/jr^ ^ f fry 2 - 5xv Wr 

- : 'f(2x h r h ■ r - 7xx h 2 4 27p)xj 

- .4 6 1 ■ ^r^-rrrv - 2xj,r • r fl 4- 7xx h 2 ) - §r(7xj, 3 - 3 x h r h 2 ) 

+ f(-rh 2 rk ■ r 4- lxr h 2 x h 4- 7x/ 4 2 r ■ r /t - 21xx& 3 )x] 

3M3 [( 2 xXgm ~r yyem 4" ^-4 eTn )x -r (yx em t xy ern )y 

T ~ .£2 e 77 l )zj /( 7 R) 

- /i 3 [-lox em Xhr 4- 3r h ■ r cm r - 15xr sm x„ 

4- 3r • TkTem + 3(r ■ rem)(-5x*x 4- r/i) 4- 105xXemX;,X 
- 15xx em r^ - loxrft • r em x - 15x em r • r,- t xj /(' O) 0 


1.0198 x 1C Y'CrAc 
to}) (7 — 1 )Dx 4- Th 4- ri; 3 
1.0198 x 10 17 C R Aa 


m 1 


l)Dx 4- Th ii 3 1 




[(7 4- l)Dx - r h 
4 l)Dx 4- rh] 


r j 


c(b') -s(u ,f ) 

01 


c(6) 

0 


f 1 °. 

0 

s{p) c{ib f ) 

0 


0 

I 

0 

j 0 C (<P) 


0 o' 

c 


-MB) 

0 

c(0) 

j_0 s(d>) 

c(<b) 


F thr 


thrust 

m 


> ER3B 


> lunar 


} SRP t 
| SRP/, 


thrust 


where v r = ^ 4- /. 

The following example permits comparison of the relative contribution of the terms 
to the telescope motion. The initial conditions are the same as those presented earlier in 
Table 1. As in Figure 3, Fignre 4 presents the solution to four different force models selected 
from the relative motion summary equation above. Where hub position is required, it was 
obtained from separate integration of the full, unexpanded hub equations from the previous 
section, using the same force model as for the relative motion. 

Once again, the reference solution is represented in the figure by the x-axis. This refers 
to the circular restricted model case. The other models depicted in Figure 4 represent the 
addition of ellipticity (ER3B), lunar, and solar radiation (SRP) effects over 20 days. Solar 
radiation is again treated as being applied to both the telescope and to either a real or 
phantom spacecraft at the hub, with the same physical characteristics as the telescope. 


SENSITIVITY TO HUB POSITION ERROR 

One area of concern to mission planners involves the limitations of accuracy 7 in knowledge 
of the hub position. In particular, as the relative motion equations are dependent upon the 
hub position, it is important to understand how sensitive the relative motion solutions are 
to errors in the hub position. To address this issue, a detailed linear variational analysis 
was performed, clearly demonstrating that the errors in telescope position relative to the 
hub, based on knowledge of hub position to approximately 1.7 km. are sufficient])’ small 
that they may be ignored. 

Several sample tests of this behavior were conducted. For one test case, the initial 
conditions of Table 1 were used here as a nominal set of initial conditions for the hub and 
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Figure 4 

Effects of Perturbations on Relative Distance — Expanded Equations 


telescope. The integration was then performed with the hub offset from its nominal initial 
state by 1 km in each directional component. This was to simulate an initial error in hub 
position. As seen in Figure 5, the telescope position relative to the hub differs from the 
nominal case by millimeters over 40 days. 



Figure 5 

Relative Motion Error Caused by 1.7 km Initial Hub Position Error 


The simulation was also repeated with the initial hub position offset from the nominal 
state by 17.3 km (10 km in each direction). The resulting relative motion error is roughly 
the same as that of Figure 5. scaled by a factor of 10. 
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SUMMARY AND CONCLUSION 

This report details the further work describing the formation flying between spacecraft near 
the Sun-Earth L 2 libration point, beginning with the circular restricted three-body problem 
for the hub motion about L?. 

These analyses develop the elliptical restricted three-body problem from previous work 
with circular problem [1]. The earlier in-depth work with the circular restricted problem 
was used as a basis from which to address the following perturbations: 

• elliptical orbit of Earth-Moon about Sun 

• lunar gravitational effects 

• solar radiation pressure effects 

• thrusters on vehicle 

These were incorporated as additive perturbations to the circular restricted three-body 
problem with expansions of varying levels of fidelity. For use as a baseline, the equations of 
motion were first developed in their full nonlinear form; then, they were expanded through 
appropriately small contributions. The equations were implemented through their coding in 
a MATLAB simulation. One example was presented using the earlier identification of valid 
initial conditions that excite only oscillatory motion in the linear modes. The results shown 
in Figures 3 and 4 demonstrate the dominant perturbation due to the elliptical motion of 
Earth about the sun. In this example, the perturbations due to lunar and solar radiation 
pressure are negligible. Overall, the figure indicates that during the typical 5-day science 
observation of this example, these perturbations contribute less than 1 m of relative position 
difference. The proposed electric propulsion should have little trouble maintaining position. 

The expanded circular and elliptical portions of the full model through terms which 
are linear in the coordinates of the telescope position relative to the hub and no more 
than cubic in the coordinates of the hub position relative to L 2 are presented above. The 
derivation describes the magnitude of the various terms and why some were truncated. 
Additionally, the lunar gravity model was expanded and some terms truncated due to very 
small perturbations upon the motion. As before, the perturbations due to the elliptical 
motion of the Earth about the Sun are dominant. Again, lunar and solar radiation pressure 
effects are negligible. 

This work’s uniqueness stems from its complete description of the primary perturba- 
tions to the relative motion between nearby spacecraft. In the course of the analysis, the 
dominance of the perturbing effect of the elliptical motion of the Earth about the sun was 
identified and verified. Contributions due to lunar gravity and solar radiation pressure are 
nearly negligible in the chosen examples. 
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CONSTANT PARAMETERS 


Table 4 

PHYSICAL CONSTANTS [3] 

(E-M bary. = Earth-Moon barycenter) 


gravitational parameter, Sun alone 

Ml 

132,712,440.017.987 kiiU/s 2 

gravitational parameter, Earth alone 

M2 

398,600.4415 km 3 /s 2 * 

gravitational parameter, Earth-Moon 

M2 

403,503.236 km 3 /s 2 

gravitational parameter, Moon alone 

P3 

4,902.8003 km 3 /s 2 

astronomical unit 

AU 

149,597,870.691 km 

mean E-M bary. distance from Sun 


1.000001018 AU 

eccentricity of E-M bary. orbit about Sun 

e 

0.01670862 

mean motion of E-M bary. orbit about Sun 

n 

0.199106385 xl0~ 6 rad/s 

L 2 distance ratio 

7 

0.01007824 


Table 5 

COMPUTED VALUES AND COEFFICIENTS 


reference Sun-Earth distance 

D 

149,618,905.218739 km 

gravitational coefficients 

A 3 

1.16556055765939 x 10~ 3 1/day 2 

(see Equation (9)) 

a 4 

5.84525170441422 x lO -10 l/(km-day 2 ) 


As 

3.86396147244215 x 10~ 16 l/(km 2 -day 2 ) 


A 6 

2.56240418152728 x 10~ 22 l/(km 3 -day 2 ) 
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